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Abstract 

Single top production in e + e~ annihilations is searched for in data collected by 
the L3 detector at centre-of-mass energies from 189 to 209 GeV, corresponding to 
a total integrated luminosity of 634 pb _1 . Investigating hadronic and semileptonic 
top decays, no evidence of single top production at LEP is obtained and upper 
limits on the single top cross section as a function of the centre-of-mass energy are 
derived. Limits on possible anomalous couplings, as well as on the scale of contact 
interactions responsible for single top production are determined. 
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Introduction 



One of the fundamental features of the Standard Model is that neutral currents are flavour 
diagonal, therefore any flavour changing neutral current (FCNC) process may occur only at 
second or higher orders through loops. The FCNC interactions of the top quark offer an ideal 
place to search for new physics and can be studied either in t — > qV (q = u, c; V = 7, Z) decays 
at the Tevatron or in single top production at LEP. In the Standard Model, the predicted 
rates of such processes are very small [1., but new physics beyond the Standard Model could 
lead to a significant increase. Enhancements of some orders of magnitude are predicted in 
two-Higgs-doublet models and supersymmetric models [2]. Flavour changing multiple-Higgs- 
doublet models further enhance the rates up to 10~ 5 [Sj, and models with FCNC coupled singlet 
quarks, compositeness or dynamical electroweak symmetry breaking 4] can yield a branching 
ratio of about 10~ 2 , in the reach of present colliders. The same considerations apply for the 
predicted single top production rates Experimental upper limits at the 95% confidence 
level on the FCNC branching fractions of the top quark were set by the CDF Collaboration as 
BR(t —> c(u)7) < 3.2% and BR(t — ► c(u)Z) < 33% 0. Other studies of single top production 
were carried out at LEP 

In this Letter, we consider two theoretical models for single top production e + e~ — > tc 
an interpretation involving [tce + e~] contact interactions and an approach which describes the 
process in terms of anomalous couplings. 



Single Top Production through Contact Interactions 

An effective flavour changing vertex can be parameterized through [tce + e~] contact interac- 
tions 8J. The total unpolarised cross sections for e + e~ — > tc and charge conjugate can be 
written as: 
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where Vij, Sy, Tij represent the coupling constants of vector, scalar and tensor fields, respec- 
tively (L and R stand for left-handed and right-handed fields), C = s/A 4 x (3 2 /Att(1 + (3) 3 , 
(3 — (s — m 2 )/(s + m 2 ) and A is the energy scale parameter for the process. The predicted 
single top production cross sections are shown in Figure ^1 for three different assumptions on 
the Lorentz structure of the [tce + e~] vertex. In this Letter, we give lower limits on the energy 
scale A. No directly comparable previous limits exist. 



Anomalous Couplings 

The Born-level cross section for the single top production process e + e~ — > tc in the presence of 
anomalous couplings tVc, with V = 7, Z jHHHIl, is given by: 

in the limit of a massless c quark; k 7 and k-i define the strength of the anomalous coupling for 
the current with a photon and a Z boson, respectively, s is the centre-of-mass energy squared, 
a is the fine structure constant, e q = 2/3 and m t are the charge fraction and the mass of 

Charge conjugate is assumed throughout this Letter. 
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the top quark, aw = 1 — 4sin 2 6 l w, <Sw = sin 2 2#w, #w is the Weinberg angle, M. t = m\j s 
and Aiz = m\/s. QCD corrections are taken into account following the prescriptions of 
Reference ITT1 Moreover, the effect of initial state radiation (ISR) must be included. Born-level 
and QCD+ISR corrected predictions for the cross section are shown in Figure ^p. The FCNC 
branching fraction limits set by CDF correspond to upper limits on the anomalous couplings 
of k 2 < 0.176 and k\ < 0.533 0, at the 95% confidence level. 

Data and Monte Carlo Samples 

This study is based on 634 pb _1 of data collected by the L3 detector [T2j at LEP at y/s = 
189 — 209 GeV. This integrated luminosity corresponds to the seven ranges of average centre- 
of-mass energies shown in Table ^ 



< v/i > (GeV) 


188.6 191.6 195.5 199.6 201.8 204.8 206.6 


Luminosity (pb _1 ) 


176.8 29.8 84.1 84.0 37.7 86.0 135.5 



Table 1: Average centre-of-mass energies and their corresponding integrated lumi- 
nosities. 



Crucial to this analysis is the prediction of Standard Model backgrounds, which rely on the 
following Monte Carlo (MC) programs: PYTHIA p3] and KK2f [Hj for e+e" -> qqft), KO- 
RALW [IH] for e+e" -> W+W", KORALZ [Ej for e+e" -> r+r", PHOJET [H| for two-photon 
interactions and EXCALIBUR ^B] for other four-fermion final states. Single top production 
MC events were generated with a modified version of PYTHIA [T^]. QCD colour reconnection 
effects are taken into account in the framework of the Lund fragmentation model [20J, forcing 
the top decay t — > W + b before the fragmentation takes place. 

The response of the L3 detector is simulated using the GEANT [21] program, taking into 
account the effects of multiple scattering, energy loss and showering in the detector. Hadronic 
interactions in the detector are modeled using the GHEISHA [21] program. Time dependent 
detector inefficiencies, as monitored during the data taking period, are also simulated. 



Analysis Procedures 

In the reaction e + e _ — ► tc at LEP centre-of-mass energies, the top quark is produced almost 
at rest and quickly decays via t — > W + b without forming top-flavoured hadrons. Depending 
on the subsequent decay of the W boson, the expected final state signatures are two jets, one 
lepton and missing energy (tc — > W + bc — * \ + uhc) or four jets (tc — > W + bc — > qq'bc), hereafter 
referred to as the leptonic and hadronic channels, respectively. In both channels, the c-quark 
energy E c has a fixed value for a given yfs: 

E, = ,3) 

whereas the b-quark energy has an almost fixed value which does not depend on y/s, in the 
limit of a top quark at rest in the centre-of-mass frame: 

E> - (4) 
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where m\y is the mass of the W boson. After hadronisation, the c and b quarks yield jets with 
almost fixed energies. Since both channels are characterised by the presence of a b quark, a 
clear signature exists and is exploited by using a b-jet tagging algorithm [22], which is mainly 
based on lifetime information. 

Standard search procedures are applied to both the leptonic and hadronic channels. First, 
a preselection is applied which significantly reduces the background while keeping a high signal 
efficiency; this is especially effective against background from low multiplicity events and from 
two-photon interactions. Then, a further set of channel-specific selection criteria is chosen to 
increase the signal-to-background ratio. Finally, a discriminating variable is built using a neural 
network technique. 

Leptonic Channel 

The signature for the leptonic channel is one energetic lepton, large missing momentum and 
two jets with a large difference in energy. The most energetic jet is assumed to stem from the 
hadronisation of the b quark. The e + e~ — ► W + W~ and e + e~ — > qq(7) processes constitute the 
main background. 

A preselection is applied requiring events to have at least three tracks, more than 15 calori- 
metric clusters and a visible energy greater than 0.25 x t/s, but less than 0.9 x y/s. The presence 
of a well identified lepton is required; if there is more than one reconstructed lepton, the most 
energetic one is retained. An electron candidate is identified as a track with an associated 
cluster in the electromagnetic calorimeter; a muon candidate is reconstructed as a track in 
the central tracker matched to one in the muon spectrometer; a tau lepton is identified as a 
low-multiplicity jet. All clusters in the event, except the ones associated with the reconstructed 
lepton, are combined to form two hadronic jets using the DURHAM algorithm (2H|. The jet 
axes and the missing momentum vector must be at least 15° and 26°, respectively, away from 
the beam axis. 

To further reject background events, the energy of the lepton is required to be at least 
10 GeV. The energy of the most energetic jet is required to exceed 60 GeV, whereas an upper 
bound is set on the energy E\ of the least energetic jet, as detailed in Tabled The width of the 
least energetic jet must be less than 0.4, the jet width being defined as the scalar sum of the 
transverse momenta of the jet clusters, normalised to the jet energy. The missing momentum is 
required to exceed 25 GeV, the lepton plus missing momentum invariant mass must be larger 
than 20 GeV and the two-jet invariant mass M qq - has to lie outside a window around the W 
mass, as given in Tabled The distributions of some relevant variables are shown in Figure El 
From the overall selection, 346 events are left in the data sample, with 357.0 ± 1.8 expected 
from Standard Model backgrounds. The signal efficiency is 10.6% as detailed in Table 01 which 
also gives results at all centre-of-mass energies. In order to enhance the separation between 
signal and background, a neural network technique [21] is then employed. The most important 
event variables are used as inputs to a neural network with 10 input nodes and two output 



(GeV) 


188.6 191.6 195.5 199.6 201.8 204.8 206.6 


Ei (GeV) < 
|M q q-m w | (GeV) > 


17.0 21.0 23.0 27.0 30.0 34.0 34.0 
14.0 11.0 8.5 6.0 6.0 3.0 3.0 



Table 2: Centre-of-mass energy dependent cuts used for the leptonic channel. 
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nodes, O tc , Oback, corresponding to the signal and the background, respectively. The input 
variables are related to the magnitude and direction of the missing momentum vector, the fa- 
tag value, the invariant mass of the two-jets system and the invariant mass of the lepton plus 
missing-momentum system. The final discriminating variable is then obtained as the product 
= O tc x(l-O back ). 





Leptonic channel 


Hadronic channel 


y/s (GeV) 


Data 


T ) 1 

Back. 


big. ett. [%) 


Data 


Back. 


big. ett. (%) 


1 oo a 

188.6 


10 


13.6 ± 0.3 


10.7 ± 0.2 


95 


76.6 ± 0.8 


21.0 ± 0.4 


191.6 


3 


4.6 ± 0.2 


11.1 ± 0.3 


14 


14.3 ± 0.5 


23.1 ± 0.6 


195.5 


23 


21.1 ± 0.4 


10.4 ± 0.3 


43 


38.7 ± 1.0 


22.1 ± 0.6 


199.6 


35 


40.5 ± 0.4 


10.7 ± 0.3 


37 


38.1 ± 0.5 


21.5 ± 0.6 


201.8 


22 


23.7 ± 0.2 


10.2 ± 0.3 


19 


19.4 ± 0.4 


21.6 ± 0.6 


204.8 


104 


99.1 ± 1.1 


10.5 ± 0.3 


50 


41.1 ± 0.7 


20.7 ± 0.6 


206.6 


149 


154.4 ± 0.9 


10.5 ± 0.3 


63 


59.7 ± 0.6 


20.2 ± 0.5 


All 


346 


357.0 ± 1.8 


10.6 ± 0.1 


321 


287.9 ± 1.6 


21.1 ± 0.2 



Table 3: The number of data events for each centre-of-mass energy range, the ex- 
pected number of background events and the Monte Carlo predicted signal efficien- 
cies, for both the leptonic and hadronic channels. The errors are statistical only. 
The signal efficiency is for a top quark mass of 175 GeV. The corresponding W decay 
branching fraction is taken into account in the signal efficiency. 



Hadronic Channel 

The signature for the hadronic channel is four hadronic jets, two jets having an almost fixed 
energy and two being the decay products of a W. In addition, one jet must have a strong b-jet 
signature. The e + e~ — > W + W~ and e + e~ — > qq(7) processes constitute the main background. 

A hadronic preselection, which requires a visible energy of at least 0.7 x ^/s and an effective 
centre-of-mass energy larger than 0.85 x ^/s, is applied. The effective centre-of-mass energy 
is computed after having removed photon radiation in the initial state. We also require more 
than 20 reconstructed tracks and more than two jets, built using the DURHAM algorithm with 
a resolution parameter y cut = 0.02. 

The events are forced into a four-jet topology by changing the value of y cut . The c-jet 
candidate is defined as the jet whose energy is closest to the value expected from Equation Q 
for m t = 174.3 GeV. Using the remaining three jets, the W boson is identified as the jet pair 
with invariant mass closest to the nominal W mass. The remaining jet is assumed to be the b 
jet. 

By comparing the direction of the leading c quark in MC events with the c-jet defined above, 
it is observed that the c selection purities are about 49%, 53% and 60% for ^Js = 186.6, 191.6 
and 195.5 GeV, respectively, and approach the limit of 62% for ^Js > 201.8 GeV. 

A neural network with 24 input nodes and three output nodes is then used. The neural 
network input variables are related to such jet characteristics as their energies and masses, the 
b-tag discriminant of the b jet candidate and various event shape variables. Some of these 
distributions are shown in Figures and Eb- One network output selects the signal and is 
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used as the final discriminating variable. The other two outputs, which tag e + e~ — ► qq(7) and 
e + e~ — > W + W~ events, are used to further reject the background by applying a cut of 0.1 
and 0.3, respectively. Their distributions are shown in Figures Ob and Ell. After all cuts, the 
final sample consists of 321 data events, to be compared with 287.9 ± 1.6 expected background 
events, and a 21.1% signal efficiency, as detailed in Table El 

Study of Systematic Uncertainties 

The search for single top production discussed in this Letter relies heavily on the comparison 
between the data and its associated Monte Carlo simulations. Uncertainties in these simulations 
give rise to three sources of systematic uncertainties whose effects on the single top cross 
section are shown in Table 4. First, the finite statistics of the Monte Carlo used for the 
signal and background simulations affect the determination of the signal efficiency and the 
level of background contamination. Secondly, the background cross sections are fixed in the 
interpretation of the observed events in terms of the single top cross section. The effects of a 
variation of the e + e~ — > qq(7), e + e~ — > W + W~ and e + e~ — > TIL cross sections of 2%, 0.5% and 
10%, respectively, are reported in Table 0] Thirdly, there can be differences between the actual 
and simulated detector performance, affecting the description of variables used as inputs to the 
neural network. In particular, we study the variables used for the lepton identification, the 
global event shape, the energy measurement and the b-tag. Their effects on the single top cross 
section are also given in Table EJ that also lists effects from the modeling of other variables. 
Finally, the systematic due to uncertainties in the measurement of the integrated luminosity is 
negligible. 



Source 


Aa/a (%) 


MC background statistics 


0.4 


MC signal statistics 


1.7 


e + e~ — > qq(7) modeling 


0.2 


e + e~ — > W + W~ modeling 


< 0.1 


e + e~ — > LL modeling 


0.4 


Lepton identification 


0.1 


Event shape 


0.9 


Energy resolution 


1.6 


b-tagging 


1.4 


Other variables 


0.5 


Luminosity 


< 0.1 


Signal angular distribution 


2.0 



Table 4: Relative systematic uncertainties affecting the single top cross section. 

In addition to these effects, uncertainties in the modeling of the signal process can affect our 
results. To quantify this, various signal samples with different final-state angular distributions 
are simulated, inside the limits allowed by the anomalous coupling scenario, with the effects 
reported in Table 4. Systematic uncertainties on the signal Monte Carlo statistic are propagated 
to the signal efficiency Uncertainties on the background Monte Carlo statistic and modeling 
are accounted for by lowering the background expectations accordingly. Finally, the uncertainty 
in the detector description is propagated by repeating the analysis with the distributions of the 
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variables entering the neural networks smeared to agree with the distributions observed in data. 
An additional source of uncertainty could be the value of m t used in the simulation. The low 
momenta available for the top system at our centre-of-mass energies would imply a change in 
the event kinematics and hence in the selection efficiency Rather than assigning a systematic 
uncertainty from this source, all results are parametrised in terms of m t . 

Results 

Leptonic and hadronic final discriminating variables are shown in Figure 0] A very good 
discrimination between signal and background is achieved. No significant deviation from the 
Standard Model background expectation is observed. Combined 95% confidence level upper 
limits on the single top total cross section are derived 25] and listed in Table El For these 
limits, the branching ratio for the top decay is assumed to be saturated by t — > Wb. The limits 
are obtained for the signal process e + e~ — > tc, a deterioration of the limits of about 10% is 
found for the corresponding process e + e~ — > tu. 

Table El lists the 95% confidence level lower limit on the energy scale parameter A of a 
possible [tce + e~] contact interaction, described in Equation (JTJ). Referring to the anomalous 
coupling formalism described in Reference [HI using the cross section expression of Equation (J2J), 
an exclusion region in the k z vs. k 7 plane is obtained, as displayed in Figure Et- QCD and ISR 
corrections as well as flavour changing decays of the top quark through anomalous vertices are 
taken into account in the limits computation. Using the Born-level cross section of Equation ((2j), 
a corresponding exclusion region in the BR(t — > Zq) vs. BR(t — > jq) plane is found, as shown 
in Figure Eb- The anomalous coupling formalism upper limits are summarised in Table 

In conclusion, no evidence for single top production at LEP is observed and possible new 
physics responsible for this process is constrained. 



< V* > 




(GeV) 


188.6 


191.6 


195.5 


199.6 


201.8 


204.8 


206.6 


0"95Ot = 


170 GeV) 


(pb) 


0.36 


0.87 


0.77 


0.54 


0.62 


0.63 


0.51 


0"95 {mt = 


175 GeV) 


(pb) 


0.22 


0.73 


0.67 


0.45 


0.56 


0.51 


0.39 


0-95 (mt = 


180 GeV) 


(pb) 


0.21 


0.75 


0.64 


0.42 


0.52 


0.48 


0.37 



Table 5: Measured 95% confidence level upper limits, cr 95 , on the total cross section 
for single top production as a function of the centre-of-mass energy. The limits are 
given for three assumptions on the top quark mass. 
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A (TeV) 






Vector coupling 


Scalar coupling 


Tensor coupling 


m t = 


ill) Ore V 


U. i u 


u.uo 


1 9zL 


m l — 


175 GpV 


0.75 


0.65 


1.24 


m t = 


180 GeV 


0.70 


0.60 


1.16 



Table 6: Measured 95% confidence level lower limits on the energy scale parameter A 
in TeV. Three different scenarios for the coupling constants are considered: vectorial 
(Vij = 1, S RR = 0, T RR = 0), scalar (Vij = 0, S RR = 1, T RR = 0) and tensorial 
(Vij = 0, S RR = 0, T RR = 1). Limits are given for three values of the top quark 
mass. 



mt (GeV) 


170 


175 


180 


|«z| 


0.38 


0.37 


0.43 


1 1 


0.43 


0.43 


0.49 


BR(t -> Zq) 


13.6% 


13.7% 


17.0% 


BR(t - 7 q) 


4.4% 


4.1% 


4.9% 



Table 7: Measured 95% confidence level upper limits on the anomalous-coupling 
parameters kz and k 1 and on the FCNC top decay branching fractions. Limits are 
given for three values of the top quark mass. 
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Figure 1: Theoretical total cross section for single top production as a function of 
the centre-of-mass energy for three different top quark mass values for (a) [tce + e~] 
contact interactions with different assumptions on the Lorentz structure and an 
energy scale parameter A = 1 TeV and (b) the model described in Equation (J2J) 
where CDF limits on the values for the anomalous constants are assumed. 
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Figure 2: Distributions for the leptonic channel of (a) the b-tag discriminant variable 
for the most energetic jet, and (b) the energy E\ of the least energetic jet. The signal 
histogram is for a top quark mass of 175 GeV and is normalised to the number of 
data events. Background expectations are also shown. 
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Figure 3: Distributions for the hadronic channel after preselection of (a) the b-tag 
variable of the b jet, (b) the event thrust and (c), (d) the neural network outputs 
related to e + e~ — > qq(7) and e + e~ — > W + W~. The signal histogram is for a top 
quark mass of 175 GeV and is normalised to the number of data events. Background 
expectations are also shown. 
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Figure 4: Distribution of the final discriminating variable for (a) the leptonic channel 
and (b) the hadronic channel. The signal histogram is for a top quark mass of 
175 GeV and is normalised to the number of data events. Background expectations 
are also shown. 
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Figure 5: Exclusion regions at the 95% confidence level in (a) the k z vs. k 7 plane 
and (b) the BR(t — > Zq) vs. BR(t — > jq) plane for three different values of the top 
quark mass. The CDF exclusion domain is also shown. On Figure (a) the curves 
m t = 170 GeV and m t = 175 GeV are almost overlapping. 
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